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** Abstract ; We consider the non-Hnear random vibration 
model described by the system of differential equations . 
The solution can be found by using the inverse Laplace transformation . 

** Subjects: Vibration Mechanics , The Differential equations . 

NOTE: 

This worksheet demonstrates Maple's capabilities in the design and finding 
the numerical solution of the non-linear vibration system . 
All rights reserved. Copying or transmitting of this material without the 
permission of the authors is not allowed . 



1. Model Definition . 

Voltage supply ; V 

Switch : S 

G eartrain : G 

M otor : M . 

Drum : D 

Cable ; C 

Elevator cage ; E 

(Fig.l) 

The cage E is hoisted by cable C wound over a drum D driven through a gear-set G by an electric motor M 

A ssume that ; 



* the cable C has low internal friction . 

* the variation of weight supported with length change may be ignored 

* weight is concentrated in tiie cage E 

* variation of cable C with length change may be ignored 

* the cable C internal damping may be neglected . 

* drum D and gear G inertia may be ignored . 

* DC electric motor with constant magnetic field . 

* motor armature resistance & inductance may be neglected 

* relay resistance may be neglected 

* internal resistance of the voltage supply may be neglected . 
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( Fig.l ) 
( Fig.2) 



2 . The system of differential equations : 



m^ ^ = ^ (x- - X£.) - m^.g 

0)^ =—^ : K^ : coefflent of motor tran sduction . 

dx- 
dt 



r ^2 



el :=m. 



de'=<« 



= k^{x^{t)-x^{t)) + m^g e2:=aj^{t) 






/2:=«d = "g«m % ==^ « 



M 



e3:=j^xi,t) = r^(sy^ 



el := m^ 



/^2 ^ 



c,^'^''' 



= k^(x^{t)-x^{t)) + jn^g e2:=cOj^(t) = -^ e3 :=^x,(t) = r^ ng(t)cOj^^(t) 
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r 

'm 
£]:=-m^(D(x^)(0) + sx^(0)) + m^s2l^ace(x^(tXts) = 



^m(^) . d 



kp laplace(>cj( t), t s)-k^ laplace(Xj(t), t, s) 



+ ■ 



s 



E2:=a>(t) = 



eAt) 



^M 



M 



£3 :=s laplace(x^(t), t s) - x^(0) = r^^ l^lace(n^(t) co^(t). t s) 



x^(0):=0 x^(0):=0 D(x^)(0):=0 



7 ^E 5 

£"] •=^£^ l^ace(x^(t). t s) = k^ laplace(xj(t), t s) -k^ laplace(x^(t), t, s) + 



£2 := laplace( co ( t), t s) = 



s 
laplace(e„(t), t, s) 



M 



£3 := s laplace( Xj( t), t, s) = r^ lapLaodi n^ ( t) cOj^^( t), t, s) 

laplace( Xj( t), t, s) r^ l^ace( n^ ( t) (Oj^( t), t s) 
laplace( e^^ ( t), t, s) ~ s laplace( gj^ ( t), t, s ) 

laplace( Xj ( t), t, s) r^ l^ace( n^it)G>^it),t,s) 



laplace( e^^l t), t s) s K^^ l^ace( Oj^ ( t ). t s ) 



eql:={m^s' + k^)x^=k^x^ ec^:=sx= eq3 ■= 
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m. 



e, :=0.05 L :=1.5 K„ :=3.2 m, := 100 r, :=2.1 n, :=4 

£ C M £ D G 



7.875000000 e„ 2.625000000 e, 

^ s(200. s2+3.) ' s 

XOg := -0.1312500000 cos( 0.1224744872 t) + 0.1312500000 

Xg := -0.1312500000 cos( 0.1224744872 t) + 0.1312500000 



Xg := -0.1312500000 cos( 0.1224744872 lO'') + 0.1312500000 

X^ := 0.1312500000 XOl^ := 0.01607477644 sin( 0.1224744872 t) 
XL := 0.01607477644 sin( 0.1224744872 10* ) ho:= 2.65 
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The animate graphic of X and X 1 
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ylO :=t^ -0.1312500000 cos( 0.1224744872 t) + 0.1312500000 

yl := 50 ylO 



ylO := t^ 0.1312500000 - 0.1904751818 10"^ g(-oo9999985488t) cos(ioi.6658374 1) 
-0.5620631121 10' e'"""""'^*" " sii< 101.6658374 1)- 0.5571722761 10'^" (0) 

(0.50438898 10>> e' "^'^^^^^^ ^' co9(101.6658374 1) 

-0.17093024 10^*6 "•"'''''«^'«« ^' sii<101.6658374t)) 7- 0.5571722761 lO'^" (0) 



(-0.50438898 10^^ e 



+ 0.17093024 10'* e 



-0.09999985488 1) 



-0.09999985488 



cos( 101.6658374 1) 
sii< 101.6658374 1)) 7 



0.1312498095 e 



(-0.2902490119 10"" t) 



yl ■= 500 ]/10 



>mohinh:=proc(M) 
> mohinh(2) ; 



■ucF-do.M-do.b-wang.c 1 -xam.cS-xanh.TRAWIONGCO 




X^:= 



^-^"^ 



e r 



d"g h 



^ Rr^ n„ k„ s RC k„ s k„ 

iC.. ' ^E ^E 



V 



M 



/ 



e^ := 0.05 k^ := 1.5 K^ := 3.2 m^ := 100 r^ := 2.1 n^ :=4 e := 0.05 RC s + 0.05 

0.03937500000 ( 0.05 RCs + 0.05 ) 



^£==- 



.3 , r,2 



s{RCs' + s'+ 10.33593750 R s + 0.01500000000 RCs + 0.01500000000 ) 



XO^ :=0.1312500000 - 3.281250000 

E 



(_at) 



y _a = RootOf( 48 + 3200 RC _z'^ + 3200 _Z^ + ( 33075 R + 48 R C ) _Z ) 



e - ( 128 _a + K ( 1323 + 128 _a^ C ) ) 
9600 i?C a^+6400 a+ 33075 K + 48i? C 



R := 5000 



1000 



wo := 2.5 
Warning, the name changecoords has been redefined 
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t= _S7143e-1 
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Disclaimer: While every effort has been made to validate the solutions in 
this worksheet the author is not responsible for any errors contained and 
are not liable for any damages resulting from the use of this material. 
Legal Notice: The copyright for this application is owned by the author . 
The application is intended to demonstrate the use of Maple to solve a 
particular problem. It has been made available for product evaluation 
purposes only and may not be used in any other context without the express 
permission of the author . 
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